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Abstract. This paper explores a class of non-linear constitutive relations for 
materials with memory in the framework of covariant macroscopic Maxwell 
theory. Based on earlier models for the response of hysteretic ferromagnetic 
materials to prescribed slowly varying magnetic background fields, generalized 
models are explored that are applicable to accelerating hysteretic magneto-electric 
substances coupled self-consistently to Maxwell fields. Using a parameterized 
model consistent with experimental data for a particular material that exhibits 
purely ferroelectric hysteresis when at rest in a slowly varying electric field, a 
constitutive model is constructed that permits a numerical analysis of its response 
to a driven harmonic electromagnetic field in a rectangular cavity. This response 
is then contrasted with its predicted response when set in uniform rotary motion 
in the cavity. 
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1. Introduction 

In recent years developments in technology have enabled the synthesis of new materials 
with interesting mechanical and electromagnetic properties. These have, in turn, led 
to new advances in technology as well as challenges in our understanding of physics 
at the interface between classical and quantum behaviours. In particular many meta- 
materials on the mesoscopic scale have a rich electromagnetic phenomenology. In the 
presence of high frequency or high intensity electromagnetic fields many exhibit local 
or non-local non-linear electromagnetic constitutive properties. At the other extreme 
some materials exhibit a delayed response to slowly varying electric or magnetic 
fields. Ferromagnetic media respond with a delayed magnetisation while ferroelectrics 
respond with a delayed electrical polarisation. They also exhibit an ability to maintain 
a saturated internal magnetisation (ferromagnets) or internal electric polarisation 
(electrets). Anisotropic magneto- electric materials also exist that sustain both types 
of fields. When the internally induced fields are not uniquely determined by any 
externally applied field one often says that the material exhibits electromagnetic 
constitutive properties with memory. When the external fields vary periodically with 
time and the induced fields respond periodically, the process is often referred to as 
hysteretic and the corresponding non-linear constitutive relation between them may 
exhibit a discontinuous or branched structure to account for this [T] , |2] , [3] . However 
such terminology is often restricted to processes where the induced fields in the medium 
saturate at some level and where the shape of the resulting hysteresis loop, obtained 
by displaying the magnitude of the drive field against the magnitude of the induced 
magnetisation or polarisation is independent of time. When the drive field is aperiodic 
in time the memory effects may exhibit a more complex Lissajoux structure in the 
hysteretic response, particularly if the time dependent drive field contains more than 
one dominant harmonic component. 

A rapidly varying time-harmonic electromagnetic drive field may also induce 
both electric and magnetic polarisations with magnitudes dependent on the magneto- 
electric susceptibility 3-tensors of the medium. For materials with memory such 
susceptibilities will depend non-lincarly on the electromagnetic field in the medium 
and the resulting hysteretic response will involve both induced polarisations. In 
addition to hysteretic responses all materials exhibit spatial and temporal dispersion 
to some degree and may also sustain induced electric currents as a result of their 
conductivity. Even in homogeneous and isotropic media the detailed description of 
such media in terms of their basic constituents and micro-structure is rarely possible 
and recourse to a parametrised model becomes necessary [3]. The parameters of 
such a phenomenological model are sought from experiment over some range and the 
model tentatively extrapolated outside that range. The degree of extrapolation is 
often dictated by comparing the model with experiment. 

In this article a model is constructed that can describe a rigid non-dispersive, 
rapidly, uniformly rotating, hysteretic medium in an external time-harmonic 
electromagnetic field given its behaviour at rest. The model assigned to the 
medium at rest is motivated by a non-covariant model constructed by Coleman and 
Hodgdon [5] , [6] . By exploiting the inherent spacetime covariance of the macroscopic 
Maxwell equations such a model can be embedded in a covariant formulation and 
coupled naturally to such equations. This system can then be reduced to a coupled 
differential system in terms of electromagnetic fields and spatial tensors describing 
the magnetisation and polarisation in an arbitrarily moving medium. To illustrate 
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how such a system can have practical imphcations these equations are attacked 
numerically for an axially rotating hysteretic ferroelectric in a perfectly conducting 
cavity containing electromagnetic fields driven by an external harmonic electric 
current. 

2. Macroscopic Covariant Electrodynamics 

A theory will be said to admit a spacetime covariant formulation if it can be expressed 
in terms of tensor field equations on spacetime. The covariant theory of macroscopic 
electrodynamics benefits from a formulation in terms of differential forms. Aside from 
its elegance it makes precise a number of conceptual terms used in the interpretation 
of the theory of accelerated media and provides a suite of economical tools that 
streamline calculations. Such tools include the exterior product, exterior derivative, 
Lie derivative, covariant derivative, interior derivative and Hodge map [7]. These 
operations find their natural setting on arbitrary manifolds of arbitrary dimension. In 
this article they are employed on 3-dimensional Euclidean space and 4-diniensional 
Minkowsi spacetime. In the former case the Hodge map is denoted by # and satisfies 
the rule = 1 when acting on all (time-dependent) differential forms on space. 
The exterior derivative on such forms is denoted d and satisfies dd — 0. These 
two identities suffice to determine the many interrelations between the curl operator 
and div operator in Euclidean space. In spacetime the Hodge map 

is denoted by * and satisfies irk ~ (— 1)^+^ when acting on p— forms on spacetime. 
The exterior derivative of such p— forms on spacetime is denoted d and satisfies 
dd = 0. Although gravitation is neglegible in the following the relations between the 
frame-dependent spatial description of electromagnetism and its frame-independent 
spacetime description and the Hodge map that enters via the constitutive modelling 
require a metric tensor field g for their formulation. To this end one introduces a set 
of independent cobasis 1— forms , , e'^ , e'^ on spacetime and writes the Minkowski 
metric tensor field 

g== -e°(g)e°-)-g (1) 

3 

where g = ^ e'^ (g) e'^ is the induced metric tensor on Euclidean space. In terms 

~ k=l 

of these forms, ★I = e° A A A e^, and #1 = A A e^. The macroscopic 
Maxwell system on spacetime is defined in terms of the electromagnetic 2— form F, a 
polarisation 2— form H and a 4— current 3— form J: 



dF =0 (2) 

dirG^J (3) 

where 

G = eoF + H (4) 

in terms of the permittivity of free space eo- It is the responsibility of constitutive 
theory to provide auxiliary equations to render an augmented system deterministic. 
Thus constitutive auxiliary conditions should specify the dependence of H and the 
4— current J (subject to 4— current conservation, dj — 0), on F and its possible 
derivatives. Such conditions may be local or non-local, linear or non-linear, algebraic 
or differential. 
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An arbitrary observer in spacetime can be associated with the integral curve of 
an arbitrary future pointing unit time-like vector field V: g{V,V) = —1. Given g, 
such a vector field determines the f --form V by 

V = g{V, -) (5) 

It is convenient to refer to such a as a frame (in spacetime) since it determines 
the components of spacetime tensors measured by the associated observer. Thus the 
2— forms F,G,Il admit the orthogonal decompositions with respect to V: 

F ^e^ AV - cB^ ^e^ AV + * (cb^ A v'j (6) 
G^d"" AV-^^d"" AV + ir(^^AV^ (7) 

n=p^ AV-^^p"" AV + *(^^AV^ (8) 

J = AV- (9) 

c 

where iye^ = 0, iyd^ = 0, iyp^ = 0, iyJ^ = , iyB^ = 0, iyH^ = 0, iyM^ = 0, 
= 0- In these expressions e^, , denote a (time-dependent) spatial electric 
field 1— form and 2— form respectively; d^ , , denote a (time-dependent) spatial 
electric displacement field 1— form and 2— form respectively; , , denote a (time- 
dependent) spatial magnetic induction field 1— form and 2— form respectively; , 
, denote a (time-dependent) spatial magnetic field field 1— form and 2— form 
respectively; p^ , , denote a (time-dependent) spatial electric polarisation field 
field 1— form and 2— form respectively; mX , , denote a (time-dependent) spatial 
magnetic polarisation field field 1— form and 2— form respectively; denotes a (time- 
dependent) spatial 3— current 2— form and denotes a (time-dependent) spatial 
electric charge density 0— form. It should be stressed that these decompositions are 
defined for arbitrary observer fields, including those describing accelerated frames. 

In Minkowski spacetime there exist global coordinates t, x, y, z in which the above 
g-orthonormal cobasis takes the form: 

e° = cdt, e^ = dx, e^ = dy, = dz (10) 

The history of an inertial observer is then part of an integral curve of the vector field 

u = -dt (11) 

c 

and 

[/= -cdt (12) 

In such a spacetime coordinate system the inertial frame time rate of change of any 
time dependent p— form a is defined as the Lie derivative cLu a and denoted a. 

If one makes the orthogonal decompositions above with respect to such an inertial 
(laboratory) frame the Maxwell system ([2|) yields: 



^d#b^ = Q or V-B^ = (13) 
#de^ = - \F or V X E'^ = -B^ (14) 
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#d#d^ = or V-T>^ ^p^ (15) 

=d^+f or VxH^ = -D^+J^ (16) 

where = and a = Cgta. The bold-face characters E^, D*^, B'^, H^, J*^ 

refer to the traditional Cartesian components of the associated 3— vector fields in the 
inertial (laboratory) frame U. The relation G = cqF + 11 yields 

D'^ = eoE'^ + P^ (17) 

H^ = J_B^ + M^ (18) 
Mo 

and substituting these into the Maxwell equations (|14ll6p gives, 

^ - —Me" + m^' or = - —V x + (19) 

Mo Mo 

iJJ ^}_U^^j^u _^u _pu^ orE^ = l(vxH^-J^-P^) (20) 
Co Co ^ ^ 

3. Covariant Constitutive Models 

To set the models to be discussed in context it is worth recalling the simplest covariant 
constitutive model describing a non-dispersive, non-hysteretic, homogeneous, isotropic 
linear material with arbitrary 4— velocity V: 

G ^ eaier- —\ivF ^V +—F (21) 

V Mr / Mr 

where and /i^ are dimensionless constants. In the co-moving frame V these yield 

D^ = eoe.E^, = (22) 

MoMr 

with electric displacement in the co-moving frame proportional to electric field in that 
frame and magnetic field in the co-moving frame proportional to magnetic induction 
in that frame. If G and F in pip are decomposed with respect to different frames 
(possibly accelerating) then the relations between the electric and magnetic fields 
in different frames are different and involve the instantaneous relative 3— velocities 
between the frames. This is often referred to as a motion induced magneto-electric 
effect. _ _ 

Since H = G - epF, ★((jy G) ^V) = -{I + V ^ iv)G, and iyG = e^erivF 
equation (PT|) may be written as an algebraic local relation between 11, F, G and V . 

n = eo (tr -l)ivFhV^ {fir - 1) * [{iv * G) A y) (23) 

If the medium is conducting a constitutive relation for a conductivity current 
is required. In an inertial frame U with local coordinates t, x, y, z suppose the total 
current density 1— form is 

f{x,y,z,t) = j^^^j^{x,y,z,t) + j^^^{x,y,z,t) (24) 

where j^j(a:, y, 0, t) denotes a prescribed external current. A simple isotropic Ohmic 
conductivity current arises from the temporal non-local relation: 

j^,„^(x, y, z,t)= f Kit' - t)e^(a;, y, z, t')dt' (25) 
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since its Fourier transform with respect to t yields: 

i^„„,(^,y,z,c.^) -a(c.^)e^(a;,y,z,w^) (26) 

where the scalar conductivity (t{uj^) is the Fourier transform of the spatially 
homogeneous scalar K{t), is the Fourier transform of and is the Fourier 
transform of e'^. A spacetime model therefore requires a model for k. In practice one 
finds data for cr(w'^) over some restricted range of which is rarely sufficient, in 
general, to fully re-construct k by Fourier inversion. However in circumstances where 
the electromagnetic fields have Fourier components dominantly in the frequency range 
of relevance the approximation 

j^^^aix, y, z, t) ^ a(w^)e^(x, y, z, t) (27) 

often suffices. In this approximation (t{uj^) is regarded as a homogeneous scalar field 
in space. 

To extend these models to a spatially non-dispersive but anisotropic hysteretic, 
purely ferromagnetic or ferreolectric medium it is necessary to accommodate non- 
local differential constitutive relations with memory. Such a medium has a preferred 
material eigen-basis determined by its ferromagnetic or ferreolectric susceptibility 
tensor. This determines the direction in space of the induced magnetic or electric 
polaristion relative to the direction of an applied field. A model for a uni-dimensional, 
rate-independent ferromagnetic medium was constructed by Coleman and Hodgdon 
[S], [^. In the notation of this article, with inertial spatial field components taken in 
a preferred g-ortho-normal eigen-basis it took the form: 



{f,{h'[)~b'^)+h'^g,{h'^) for J- = 1,2,3 (28) 



where the magnetic field h^j and the magnetic induction were continuous slowly 
varying real- valued functions of time with piecewise continuous derivatives, , , 
positive constants aj and fj and gj specified real-valued functions on the real line. 
Following the approach described in the introduction one seeks a covariant extension 
applicable to ferroelectrics and ferromagnets that reduces for slowly varying electric 
or magnetic fields in any inertial frame to a similar type of Coleman-Hodgdon model. 
Since such a moving ferromagnet can acquire an induced electric polarisation and 
such a moving ferroelectric can acquire a magnetic polarisation, a more general inertial 
model can be constructed that exhibits a hysteretic intrinsic magneto-electric response 
to slowly varying fields in any inertial frame. Needless to say there can be no unique 
extension to accelerating media. 

To construct such covariant models we introduce a number of rank (1,1) 
orthogonal spatial tensors on spacetime. Each such tensor x rnaps 1— forms to 
1— forms and satisfies xi^i ~) — ^ s-iid x(~j^) = for any time-like vector field V. 
It has an associated g-ortho-normal eigenbasis {£^1,^2, ^^3} and in this basis 
3 

X^Y.^j£j(^£j (29) 

In terms of four (in general distinct) such spatial tensors an intrinsic magneto- electric 
hysteretic model that has the required inertial behaviour takes the form 

VvH = Xpe {^vivF) hV + xlh (Vviy ^G)^V~ 

* {xlh (Vy*y *G)AV + xle {^vivF) A v] (30) 
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where Vy denotes a covariant derivative with respect to V and the orthogonal 
spatial tensors Xq satisfy Xq{^j~) = and Xq(~i^) = with Q—pe,ph,nih,me. 
They characterize the hysteretic properties of the medium. If all these tensors are 
non-zero for all V then the material is totally intrinsically magneto-electric. 

Since Vyll ^ VyTk-II when the (Levi-Civita) 4— acceleration VyV of the medium 
is non-zero, there exists a physically distinct "dual" model 

Vv*n = xLh (Vyiy ^G)^V + xle (Vy^yF) ^V + 

* {Xpe (VyiyF) ^V + x'^h (Vyiy ^G)Av] (31) 

that has the same behavior for media at rest in all inertial frames. 

The hysteretic behavior of the medium is determined by parameterising the 
components of each spatial tensor in its associated material eigenbasis. Thus with 

xl-txlM^^J (32) 
the scalar components Xg j i^J ^Pj ^^J ^''^J) written 



Apej 



Xph,j 



A.me,j 



mh,j 



sgn {"UphJ 



<''sgn«,^^. 



sg 



7"sgn I i£vVv 



c 



nmh^ 



sgn i^vWy- 



,gn 



^ 1 

me. J J 



with e^ = ivF, = -ciy * G for j = 1, 2, 3 and constants ,0f £ 
The sgn function defined by 

( +1 for z>0 
sgn(z) — z/\z\ = < for z = 
[ -1 for z <0 

accommodates the branching during the hysteretic process 0, 

=eo(ar/r(^r)-^r^r) 



me, J 
mh.,j 



eo [aj fj {ej ) -f ) 



V 



V 



with constants , G K.+ and frame-dependent scalars 

icvivF, e 

7 

cigviv * G 
= igviy^, rn 



cigviv * G, 



I The factors of eo imply that ^Y^^ t/<so, ^, ^me.!/"^"' - ^''^ dimensionless. 



(33) 
(34) 

(35) 
(36) 



(37) 

(38) 

(39) 

(40) 
(41) 

(42) 
(43) 
(44) 
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The real valued function characterizes the details of the hysteretic process. 
For the explicit computation below it is taken to have the form 

ff{z) = tanh{/3fz) (45) 

with I3f e R+ and where lim^^±oo /f (2) = ±1. The parameters af,Pf,^f 
must be motivated by data in the laboratory frame. 

The use of the covariant derivative in (1301) ensures that a hysteretic process driven 
by a harmonic external field can give rise to a autonomous differential equation on 
any smooth branch of a hysteresis loop as described in the next section. In any inertial 
frame U both ([30]) and IpTj) yield the coupled system of partial differential equations: 

(46) 
(47) 



1.6. 



where |ej| = Iz^vVy 



pe - U \ 
3 1 


pej 1 ^ 


-^f ley 






1 

C 




ph,3 1 ^ ''j 




ph-J 


mhiU 
3 "j 


1*^. ■ 

1 mn,,j 






mn.j 


4 


-T'^ 1 


me, J \ 






me J 


and 




isvVv- 











(48) 
(49) 



4. A Particular Model at Rest in an Inertial (Laboratory) Frame 

We consider the special case of a non-magneto-electric medium (|30|) where %, 
and Xme = and: 



ph 



VvH = Xpe (Vl.^yF) A - * (^X™„ i^v^v *G)AV 
In terms of spatial fields it follows from ([8]) that (|50|) may be written 



(50) 



VvP^ AV + /\\/vV - 



AV 



(51) 



For the medium at rest in the inertial frame U — ^dt one has V = U and equation 
(ISTI) becomes 



p'^ A + * 

where for any spatial p— form ^ 



m 
c 

(7 



U 



AC/ =Xpe(e'') AC/ + . 




A t/ 



(52) 
(53) 

since the (laboratory) inertial cobasis is parallel. Taking components in the direction 
of U and its orthogonal subspace yields 

-xye(^^) (54) 



p 



X-mh 



(55) 
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From to one has in the U frame decoupled equations for the hysteretic electric 
and magnetic polarisation fields in the medium: 

=^^re7K,|+0f (56) 

(57) 

where p'^ — pi dx + P2dy + p^ dz etc. 

Over an arbitrary time interval each of the above differential equations is a 
non-autonomous evolution equation at each point in the medium describing induced 
polarisations as a function of time. However in certain time domains their evolution 
is controlled by autonomous ordinary differential equations. When the drive fields are 
harmonic in time the resulting solutions to such equations may exhibit a limit cycle, 
(often identified as a hysteresis loop) composed of piecewise smooth line segments. For 
such cycles controlled by ()56|) the location and number of piecewise smooth segments 
in the cycle (branches) are determined by the location and number of simultaneous 
zeroes of e^ and p^ in the cycle. Similarly in cycles controlled by ()57|) branches are 
determined by the location and number of simultaneous zeroes of and in the 
cycle. The choice of parameters in the function determines such locations and also 
the degree of induced saturation in each branch of a limit cycle during the process. 
When e^ and p^ have more general time-dependences the above equations give rise to 
solutions that may exhibit hysteretic loci containing self-intersections and/or no limit 
cycle. 

Thus at each spatial point, (j56p can be written 
6'^ dv^ 

W = ^= -r*P^e.sgn(*,^,,,.) + ersgn(e7)vl/^^,, (58) 

For each of the possible combinations of sgn(ey) and sgn(\E'pg ^), ((58)) yields a first- 
order differential equation for pf{e) with branched solutions through any point 
{Pq ,Eq). If one denotes sgn{ef) hy Se and sgn{'^^) by S^i, and furthermore assumes 
that parameters are chosen so that Sxs, is constant on a particular branch a solution 
through this point can be written 

p'^ief; P^, E^,S^, Se) = P^ e~^'^'^^ -^^^ + 

e-vTiT-f^P^aP; / ' /f (T) e"r«rTdT (59) 

where rj^^ = eo (k^'^S'* -f O^Se)- Although Se is -1-1 when (e^) is increasing with 
time or —1 when (e^) is decreasing with time the sign of S<ii will in general depend 

on the state p^{ej). Furthermore when the parameters Oj , f3^ , , Kj and 6^ vary 
inhomogeneously with position in the medium this state will also depend explicitly on 
position. In such situations an analytic solution is no longer in general possible. 
Analogous solutions may be written for ([571) and when k™^ — such branched 
solutions coincide with those in [B] when 9J^^ is constant. 

5. A Particular Model Rotating in an Inertial (Laboratory) Frame 

In this section the model is applied to a rigid uniformly axially rotating cylinder 
of radius R in an external time harmonic electromagnetic field in an effort to see 



rnh.j I 
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the order of magnitude of effects on a hysteretic process in an accelerating medium. 
The differential constitutive relations are appended to the Maxwell system and the 
resulting coupled partial differential system analyzed numerically. The system is 
wholly enclosed in a rectangular 3-dimensional computational domain that simulates 
the rotating cylinder in a perfectly conducting cavity. The choice of parameters in 
the model is motivated by establishing results that may be compared with hysteretic 
behaviour of a medium at rest in a low frequency electromagnetic environment. The 
modern measurement of ferromagnetic and ferroelectric hysteresis in such situations 
is an experimental art and owes much to modern digital and piezoelectric technology. 
Needless to say great care is exercised to accommodate effects such as sample geometry, 
material electrical conductivity variations with frequency, thermal drift instabilities 
and material spatial inhomogeneities [5]. In this article we will use a particular 
hysteresis measurement [5] more as a guide to a reasonable parameter set for our 
model rather than a targeted fit to a specific material specimen. 

The model is applied to a 1.36m radius cylinder, rotating at an angular 
speed 1497rpm inside a perfectly conducting cavity with sides of lengths 
5.45mx5.45mx2.18m. The cavity fields are driven by an assembly of externally 
prescribed currents inside the cavity at 250Hz, well below the lowest natural 
electromagnetic mode of such a cavity. For the above angular speed and cavity 
dimension a lattice 20 x 20 x 8 discretization (using a fast workstation) permits 
exploration of the hysteretic processes in the rotating medium when the above 
differential equations are coupled to the macroscopic Maxwell equations. In this 
approach all components of spatial electromagnetic fields and spatial polarisations 
are with respect to the laboratory frame U and interest is directed to how these fields 
evolve parametrically with time and hence with each other in such a rotating medium. 

We consider a purely ferroelectric medium that is magnetically inert in its rest 
frame. This is simply achieved by setting each Xmh j above to zero. The first part of 
the computation involves the determination of the material eigen-basis in the rotating 
medium. To this end it is natural to adopt a polar Minkowski cobasis rather than a 
Cartesian one. With the cylinder rotating about the laboratory z— axis with uniform 
angular speed fl radians per second the Minkowski metric in terms of a cylindrical 
polar basis e*^ — cdt, e^ = c?r, = rdcj), e^ — dz, becomes 



g = —c^dt ® dt + r'^dcj) ® d4> + dr ® dr + dz ® dz 



(60) 



where 



~ + y"^ ^ X — rcoscf), y — rsincj) 
= xdy - ydx, c^r = ^ {xdx + ydy) 

d<j) = — (xdy — ydx) , dr = — [xdx + ydy) 



(61) 
(62) 



(63) 



The bulk 4-velocity field of the cylinder is. 




(64) 



Since rVl <^ c for r ^ R, 



y ~ -St + - {xdy - yd^) 
c c 

n 

y ~ — cdt H {xdy — ydx) 



(65) 



(66) 
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- 1 . By symmetry the material basis may be written 
£^ = dr, £^ = dz (67) 



and g{V, V) 




- 1 w 














where 







= d9t + /^a^ (69) 
= - c^adt + r^pd(j) (70) 

for constants a and $. Imposing the orthogonahty conditions g{W^ ,V) = and 
giW^,W^) = yields 

/3=- (71) 

All fields can now be projected into the laboratory U frame, terms of order 
{rfl)'^ /c^ removed and re-expressed in terms of the Minkowski cobasis in t,x,y,z 
coordinates. _ _ 

Since VyU = 0, VuU = the left hand side of §0^ is 

Up^ + - [xdy - yd^)p^) AU- -UM^ 



^ {xdy - yd^) (72) 
and the right hand side of ([50)) is 

XY, (yvivF) AV = xYe (Vye^) {xdy - ydx)^ (73) 



Here the six components of ([5(1)) in the basis ([TU)) yield the following coupled PDE 
Electric Sector 



system for the functions p^(i, a;, z),Py (t, x, y, z),p]! {t, x, y, z),m^{t, x, y, z), 
my{t,x,y,z),m^{t,x,y,z) 



p'^ + n {xdy - yd,^)p^ = \ {y^xYd + ^'^xYea) 



{e^ + ilxb^, + n {xdy - yd,) e^} + 

{-Xpe.l + Xpe.2} ■ 

{eY + ^yb'^ + ^{^dy~yd.)eY} (74) 



pY + n {xdy - yd,)pY = ^ {x^x^:,,! + y\Ye,2} 



{eYp,y + ^ybY + ^ {^dy - yd.) } + 
1 r j_ y \ 

■^Xy ■[ Xpe,l + Xpe,2 / • 

{e^ + r!x6^ + r!(a;9y-y5,)e^} (75) 
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p'^ + n {xdy - yd^) p'; 



,3 {e^-Q [xb^ + yb^] + n {xdy - yd,) e^] 



Magnetic Sector 



with 



mj; + O (xdy - yd,) mj^ 
rhy +^ {xdy - ydx) rUy 
rh^ + O {xdy - yd,) = - flXpe,2 [xe, + V^^] 



(76) 

(77) 
(78) 
(79) 



Xpe,i - eo < Q!l /l + $1 



sgn eo < af/f 



^ r 



efsgn 



-ye, + xCy - Q.y ( xdy - yd,) e, +Vlx {xdy - yd,) e 



S2 



xp^ + yPy , fir u 

r 



XP^+yPy , fir 

2_ + —rn'r 

r &■ 



+ y^'i + f^r-^^^ + ^x {xdy - yd,) el + ^y {xdy - yd,) e^ 



and 



X^e,s = eo a^fr (e^ - n [xb^^ + yb^] ) - ^; 



[/ _ ^xm^ + ym^ 

Pz o 



xm^ + yniy 



^fsgn (e^ - O [2/ 6^ + X 6^^] + 17 ( ar^^ - yd,) e^) } 



Xto/j,! — 0> Xm/i,2 — 0, Xm/i,3 — 



(80) 



The presence of the discontinuous sgn functions in these equations is responsible for 
the branched structure of their solutions. 
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Figure 1. The computational domain is the space inside a perfectly conducting 
rectangular cavity with sides with lengths 5.45mx5.45mx2.18m. Figure 1 shows 
an X — y cross-section of the cavity at constant z. The computation is done on a 
20 X 20 X 8 3-dimensional lattice in the cavity. The fields in the cavity are driven 
by a spatially divergence- free time-harmonic current density that flows parallel to 
the pairs oi x — y, y — z and z — x interior faces of the cavity without penetrating 
the cylindrical ferroelectric material with radius 1.36m, centered at the origin, 
oriented with its axis of rotational symmetry along the 2— axis of the cavity. A 
time-harmonic current density is taken proportional to the spatial vector (z, z, x) 
modulated to have support that excludes the cylinder. In the figure the gradation 
of shading from dark to light around the x — y perimeter depicts, from minimum to 
maximum, the magnitude of the current density in the cross section at a particular 
instant of time. 



Before coupling these constitutive equations to the Maxwell equations a further 
simplification is adopted. Suppose the medium has a uni-axial symmetry direction 
that can be aligned along the laboratory z— axis when it is at rest in the laboratory 
frame. Since the rotary motion is transverse to this axis this implies: 



1 = 

Ape,l ^7 



Xpe,2 = 



{«fsgn(eo [a^fr (e^) - e.'pl]) 



6*^ sgn (e 



yb. 



n 



Px 

P'i 



The equations (|74])-(|80l) then simplify to 
-^{yd^~xdy)v'i 
fl{ydx -xdy)py 

n {yd, - xdy)p^ + eo {a^frie}^) - CfpH • 
{4^£3^sgn {eo [af/f^ (^3'') " ^^Pa^] ) + 0^ 1^3^ | } 

{ydx - xdy)m^ + 
+ xneo4' Kfne"^) - e^'p^le"^ + 

{ydx - xdy)m^ + 



(81) 

(82) 

(83) 
(84) 

(85) 
(86) 
(87) 
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where 



xdy) mi 



ym. 



y J 



yb: 



n 



yb: 



V J 



+ Vl(xdy'- ydx)e\ 



(89) 
(90) 
(91) 



and terms of order (rfJ/c)^ consistently dropped. The time derivative of the magnetic 
field in (PTjl is found using the time derivative of (IT5|) . From ([55]) this implies that 



it is sufficient to replace with in 



it is sufficient to replace 
VLpX = Vtp^ and neX = Oe^. 



with 



and dHZ]). Similarly implies that 
, dMl) and dHZl) and, to order (rVL/cf, 




Figure 2. For small amplitude drive currents in the cavity hysteretic 

polarization loci will exhibit different features at different locations in the 
ferroelectric medium. For sufficient large amplitude drive currents all such loci 
can be driven to limit cycles that exhibit saturation. In the computations the 
model parameters do not vary with position in the medium so the resulting limit 
cycles are independent of position. A typical ferroelectric hysteresis loop on the 
left displays the evolution of the 2;— component of electrical polarization inside the 
cylinder, at rest at 2; = 1.36m. This history is induced by the z— component of 
the time-dependent electric field driven by the external initially zero 250Hz drive 
current in the cavity. The evolution proceeds from the origin where and are 
zero to a time where a limit cycle is in evidence. The parameters of the model 
(ar = 3.6 X 10*, /3r = 2.0 x 10-6m/V,?r = 1-3 x l0^m^/C,Kl^ = 0.5, = 
0.5, cr = 2.6 X 10~^S/m) have been chosen so that experimental data points (dark 
circles) from reference [9] lie close to this particular limit cycle. The points on 
the left result when the computation is repeated but with the cylinder rotating at 
1497rpm. In this simulation the changes in the pz vs hysteresis loops induced 
by rotation are imperceptible on the scales employed in the figure. 



6. Numerical Analysis 

The above constitutive equations for components of the induced polarisation in the 
rotating ferroelectric are coupled to the macroscopic Maxwell equations described in 
section 2. If we assume that initially = 0, = and = it follows from ([55)1 . 
([86|) and (|87| that the complete differential system becomes 
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and 



=n (yd, - xdy)p^ + eo {af tanh(/3f'=ej') - e^'p^ ■ 

{Kl^£^sgn{eo [al^t^nh{l3l^e^)-es'p^])+er (92) 

{ydx - xdy) + 
xneonf K tanh(/3re,^) - ifp^'M''^ + 

xOeoC [artanh(/3re,^) -a^pf] |e,^| (93) 
Q. {ydx - xdy ) + 
yf)eo«r|«rtanh(/3re,^)-a^p^ 



eoe^ 
eoe' 
eoe^ 



(V X H 



(7 _ ■(/ 

'^^x Jext,x 



(V X H^)j, - ae^ - j,^,,^^ 



J ext,z 



— (VxE^), 
Mo 

— (V X E^),, 
Mo 

— (VxE^), 
Mo 



(94) 

(95) 
(96) 
(97) 

(98) 
(99) 
(100) 



X 10" 



X 10 




2 
e^(V/m) ^^qIo 



5 
h,(A/m) ^^qB 



Figure 3. These graphs display histories of computed field components 
located at the cavity point (x, y, z) = (2.45m, 2.45m, 1.36m) (which is inside the 
rotating ferroelectric cylinder) driven by a steady external 250Hz current outside 
the cylinder. During the displayed history the cylinder performs 1 complete 
revolution. The graph on the left displays the evolution of the z— components 
6z and e^, located at this cavity point, projected onto a region of the b^-e^ plane. 
On the right a non-zero magnetization mx induced by the time-dependent hx field 
is displayed. Similar behaviour is computed for function of the magnetic 

field and indicates the presence of a rotation induced hysteretic magnetisation in 
the ferroelectric medium that increases in magnitude writh the rotation speed f2. 
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To analyze these equations numerically Yea's algorithm [Tn],[II] for standard 
materials on a standard staggered E — H mesh has been extended to include non-linear 
hysteretic media with branched solutions. Since V • B*^ and V • in the cavity are 
chosen initially zero and the external current is constructed to be divergenceless the 
full set of material Maxwell equations are accommodated. The results of this analysis 
are summarized in Figures 1, 2 and 3. 

7. Conclusions 

A class of non-linear constitutive relations for materials with memory has been 
discussed in the framework of covariant macroscopic Maxwell theory. The general 
approach enables models to be formulated for arbitrarily moving media including 
those that exhibit hysteretic responses to time varying electromagnetic fields. Using 
a particular parameterized model, consistent with experimental data for a particular 
material that exhibits purely ferroelectric hysteresis when at rest in a slowly varying 
electric field, a numerical analysis of its response to a driven harmonic electromagnetic 
field in a rectangular cavity has been performed when in different states of rotation. 
The results indicate that such a model offers a means to compute numerically the 
significance of induced hysteretic magnetisation in a ferroelectric medium as a function 
of its rotation speed and the frequency of a self-consistent electromagnetic field. 
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